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Abstract 

We use Feynman diagrams to prove a formula for the Jones polynomial of a 
link derived recently by N. Reshetikhin. This formula presents the colored Jones 
polynomial as an integral over the coadjoint orbits corresponding to the repre- 
sentations assigned to the link components. The large k limit of the integral can 
be calculated with the help of the stationary phase approximation. The Feyn- 
man rules allow us to express the phase in terms of integrals over the manifold 
and the link components. Its stationary points correspond to flat connections 
in the link complement. We conjecture a relation between the dominant part 
of the phase and Milnor's linking numbers. We check it explicitly for the triple 
and quartic numbers by comparing their expression through the Massey product 
with Feynman diagram integrals. 



1 Work supported by the National Science Foundation under Grant No. PHY-92 09978. 



1 Introduction 



Let C be an n-component link in a 3-dimensional manifold M. E. Witten presented in [1] 
the Jones polynomial of £ as a path integral over the gauge equivalence classes of SU(2) 
connection A p on M: 

Z ai _ an (M, C;k)=J [VA,\ exp (^Scs) ft Tr a Texp ^jf A^dx^j , (1.1) 
here Scs is the Chern-Simons action 

S cs = ^Tre^ / d 3 x {A,d u A p - 2 -A p A u A p ), (1.2) 
Z J M 6 

h is a "Planck's constant": 

27T 

ft = — , A;GZ, (1.3) 
k 

the trace Tr in eq. (1.2) is taken in the fundamental (2-dimensional) representation and 
Tr^Pexp (§ c . A^dx^ are the traces of holonomies along the link components Cj taken in 
the aj-dimensional representations. 

The path integral (1.1) can be calculated in the stationary phase approximation in the 
limit of large k. The stationary points of the Chern-Simons action (1.2) are flat connections 
and Witten's invariant is presented as a sum over connected pieces M. c of their moduli space 
M: 

Z au ..., an (M,C;k) = £^.UM,£;A;), 

ZitLiM, C; k) = exp L (s® + £ S^ , (1.4) 

here Scs is a Chern-Simons action of flat connections of Ai c and S^ are the quantum n-loop 
corrections to the contribution of M c . 

Suppose that M is a rational homology sphere (RHS ). Then the trivial connection is an 
isolated point in the moduli space of flat connections. Let C have only one component, so 
that it is a knot /C. In our previous paper [2] we gave a "path integral" proof of the following 
conjecture which P. Melvin and H. Morton formulated in [4] for the case of M — S 3 : 
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Proposition 1.1 The trivial connection contribution to the Jones polynomial of a knot fC 
in a RHS M can be expressed as 



Z^\M, /C; k) = Z {tT \M- k) exp 



m t 2 ^ 



aJ(a,K), (1.5) 



here v is a self-linking number ofK and J(a, K) is a function that has the following expansion 
in K~ x series: 

oo n 

J{a, K) = E D m , n a m K~ n . (1.6) 

n=0 m=0 

The dominant part of this expansion is related to the Alexander polynomial of fC: 

00 sin ( J 

vraE^a" = [ordff^M, Z)] ^— (1.7) 

the integer numbers m 2 and d are defined in [2], m 2 = d = 1 if M = S 3 . 

The virtue of eq. (1.5) is that it assembles the dominant part of the 1/K expansion of the 
Jones polynomial Z^ T \M, /C; k) into the exponential and puts a restriction on the power of a 
in the preexponential factor. This enabled us in [2] to use a stationary phase approximation 
in the Witten-Reshetikhin-Turaev surgery formula in order to derive a knot surgery formula 
for the loop corrections S^ r \ A similar formula for the Jones polynomial of a link is required 
in order to generalize this result to link surgery. However the arguments of [2] which led 
to eq. (1.5) can not be extended directly to links, because there may exist irreducible flat 
connections in the link complement with arbitrarily small holonomies along the meridians 
of link components. 

A generalization of eq. (1.5) for links was derived recently by N. Reshetikhini [5]. He 
observed that if the dimensions on in eq. (1.1) are big enough, then the representation spaces 
can be treated classically: the matrix elements of Lie algebra generators in <x,-dimensional 
representation can be substituted by functions on the coadjoint orbit of radius aj and a trace 
over the representation can be substituted by an integral over that orbit. 



2 I am indebted to N. Rcshctikhin for sharing the results of his unpublished research. 
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Proposition 1.2 Let C be an n-component link in a RHS M. Then the trivial connection 
contribution to its Jones polynomial can be expressed as a multiple integral over the SU(2) 
coadjoint orbits: 



„ n 



K d 2 3j 



mK 



^\=lkj=x V 47r \ a i\ 



exp 



x 



XI ^m(ai, . . . ,a„) (li 



m=2 



1+ E K- m P m>l {a lt 



l,m=0 



here dj are 3-dimensional vectors with fixed length 



1.9) 



and L m (di, . . . , a n ), Pj )m (oi, . . . , a n ) are homogeneous invariant (under SO(3) rotations) 



polynomials of degree M. In particular, 

n 

L 2 (ai, . . . , a n ) = E hj a i ■ a j) 
zs tae linking number of the link components Li and Cj . 



1.10) 



An example of this formula for a torus link is derived in Appendix 1 of [3]. 

In Section ^ we find a set of Feynman rules to calculate the trivial connection contribution 
to the Jones polynomial Z^\ ^ a (M, £, k). This enables us to prove the Proposition |L2| and 
to derive a property (Proposition |2l3| ) of the polynomials L m (a\, . . . ,a n ) that we will use 
in [3] in order to relate the r.h.s. of eq. (1.8) to the multivariable Alexander polynomial. In 
Section [5] we formulate a conjecture that the coefficients of the polynomials L m (Si, . . . , a n ) 
are related to Milnor's linking numbers! in a way that generalizes eq. (1.10), and present 
some circumstantial evidence that supports it. In Section f| we study the properties of flat 
connections in the link complement in order to further support our conjecture. In Sections || 
and |] we use Feynman diagrams to calculate the coefficients of the polynomials L m for 
m = 3,4. We demonstrate explicitly that they are indeed proportional to Milnor's triple 
and quartic linking numbers. 



3 I am thankful to A. Vaintrob and O. Viro for teaching me about these linking numbers. 
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We will use the following notations throughout this paper: an element v of the Lie algebra 
su(2) can be presented as v — i<J a v a , here a a are 2x2 Pauli matrices which form an orthogonal 
basis in the fundamental representation of SU(2): [<J a ,<7b] = 22e abc cr c , Tia a a b = 25 a b- Three 
components v a , 1 < a < 3 form a 3-dimensional vector v. Depending on our needs, we will 
use v , v a or v to denote the same object. 



2 Feynman Diagrams 

We are going to prove Reshetikhin's formula (1.8) for the Jones polynomial of a link by 
using Feynman diagrams in order to calculate the path integral (1.1). An introduction 
into Feynman diagram technics can be found in any textbook on quantum field theory. 
Feynman rules of the Chern-Simons theory (1.2) are described, e.g. in [6]. We also especially 
recommend [7] for a nice simple introduction into this subject. 

Our goal is to "take a logarithm" of the sum of all Feynman diagrams contributing to 
a ( M > A k). This would be easy if the link £ had no components: the logarithm would 
be equal to the sum of all connected Feynman diagrams. Therefore we have to find the 
analog of "connectivity" for the diagrams with the end-points on the link components. This 
will be the connectivity in quantum theory on the coadjoint orbit which describes the atj- 
dimensional representation of SU(2) assigned to a link component Cj. The key to defining 
the connectivity is the Campbell-Hausdorf formula^ which shows how to take a logarithm of 
the product of two noncommuting exponentials. 

Let K, be a knot in a 3-dimensional manifold M. Let A(t) be a restriction of an SU(2) 



connection onto /C: 



Ait) = A M (s(t))— , (2.1) 



4 I am indebted to N. Reshetikhin for pointing to the relevance of the Campbell-Hausdorf formula for his 
derivation of eq. (1.8). 



4 



here < t < 1 is a parametrization of /C. The trace 

Tr a Pexp (J> A^dx^ = Tr a Pexp ^ jf * A(t)dtj (2.2) 

can be considered as a partition function of the quantum theory on the coadjoint orbit of 
577(2), A(t) playing the role of an external classical field. We will try to put the trace (2.2) 
in an exponential form similar to that of eq. (1.4). 

We start with a simple particular case of the Campbell-Hausdorf formula: 

oo r> 



e v e w = exp 



(2.3) 



n=0 U - 

here v and w are elements of a Lie algebra (say S77(2)), &d v w = [v,w] and B n are Bernoulli 
numbers. 

According to the Campbell-Hausdorf formula, 

e^...e""=exp[C(^,...,^)], (2.4) 

here C(vi, . . . ,v n ) is a Lie algebra valued infinite polynomial in commutators of t>j. We 
denote as C n (vi, . . . , v n ) a part of C(vi, ... , v n ) containing the nth order monomials which 
are linear in all vf. 

C a n (v 1 ,...,v n )=C: i _ an v?...v a n n - (2-5) 
A simple corollary of eq. (2.3) is the following 

Proposition 2.1 The coefficient an is a sum of tensors coming from the diagrams like 
those of Figs. 1-4 according to the rule: a dashed line with m vertices which is depicted in 
Fig. 1, goes from the upper right to the lower left and produces a factor 

^~ 1 ^ m ~m\ ^ ^,a s(1) ,...,a s(m y ( 2 - 6 ) 

here S m is a group of permutations of numbers 1, . . . , m and the tensor /£ 01 flm is a sum of 
products of structure constants of the group 577(2): 

fb,ai,...,a m ~ (~ 2 ) e aibb 1 e a 2 b 1 b 2 ' ' ' e a m b m -ia, (2-7) 



the sum over repeated indices is, of course, implied in eq. (2.7). A contribution of the whole 
diagram is a sum of the factors (2.6) of its dashed lines over the intermediate indices. 

More generally, if we pick m elements v^, . . . ,Vi m , 1 < i\ < . . . i m < n, then the mth 
order homogeneous part of C(v\, . . . , v n ) which is linear in them, is equal to 

CI a . v^---vT m . (2.8) 

Examples 

The diagram for n = 2 is drawn in Fig. 2: 

i°aC a aua2 v?v? = ±iv u v 2 }. (2.9) 
The diagrams for n = 3 are drawn in Fig. 3: 

^aC^V^Vr = \[KV 2 ],V 3 ] + 1 (h, [v 2 ,V 3 \] + [V 2 , [ Vl ,V 3 ]]) 

= ~{[vi,v 2 ),v 3 ) + [vi,[v 2 ,v 3 })). (2.10) 
The diagrams for n = 4 are drawn in Fig. 4: 

™aC a aua2 ,a 3 , ai = K [^3,^]]] + [v h [v 3 , [v 2 ,V 4 ]\] + [v 2 , [t»i, [v 3 ,V 4 ]\] + [v 2 , [v 3 , [v^vA]} 

+ [V 3 , [Vx, [v 2 ,V 4 ]}} + [f 3 , [V 2 , [«i,U4]]]) 
+^t([Kj ^2], [V3, U4]] + K [K, U2], U4]]) 
+ 7^([K, [«2, U4]] + [V2, [[Vl, V 3 ], V 4 ]]) 

+ ^([[^2,^3], [vi,v 4 ]] + [v 1} [[v 2 ,v 3 },v 4 ]]) 

+~[[[v 1 ,v 2 ],v 3 },v 4 ). (2.11) 



The Proposition [2.1| allows us to "take a logarithm" of the parallel transport operator 



Pexp (/q 1 A(t)dtj. We introduce an "iterated" integral (see also [12]): 

/ dt 1 ...dt n {A(t 1 ),...,A(t n )}= [ A^.-.A^n). (2.12) 

JO J0<t„<...<ti<l 

Note that the iterated integral depends on a choice of the zero point of t parametrization of 
JC. 
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Proposition 2.2 The holonomy operator Pexp ( f Q l A(t)dt\ can be presented as an exponen- 
tial of an infinite sum of iterated integrals: 



Pexp Qf 1 A(t)dt) = exp ia a fj C a ai _ an £ dt, . . . dt n {A a ' {h), . . . , A a "(t n )} 



(2-13) 



Let us first present a simple "physical" proof of eq. (2.10). We split the interval < t < 1 
into many small intervals Aij with middle points £j so that 

Pexp Qf* A(t)dt) = e Mti)At le A(t 2) At 2 . . . e A(t„)At n _ (214) 

Then we apply the Campbell-Hausdorf formula to the r.h.s. of eq. (2.14) retaining only 
those terms of C(A(ti)At\, . . . , A(t n )At n ) which are at most linear in any particular A(t%). 
According to the Proposition [2~T| , such terms are given by eq. (2.8) with t> im substituted by 
A(t im ). It is easy to see that a sum of all polynomials of a given order converges to the 
iterated integral of eq. (2.13). 

To prove eq. (2.13) more rigorously we may use the following presentation of the holonomy 
operator: 

Pexpf ( l A{t)dt) = j l dt x ...dt n {A{t x ),...,A{t n )} 

00 1 r l 

E ~i E / dt i,---, dt n {A(t a (i), A(t s{n) )} 



E-r/ dt 1 ,...,dt n P[A{t l ,...,A{t n % (2.15) 
_ n n\ Jo 



n=0 •"■ seS n ' 

def. 

n=0 

here P[A(ti, . . . , A(t n )] is a path-ordered product, i.e. the Lie algebra valued forms A(ti) are 
multiplied in the order of the values of their arguments tf. for the largest U, A(U) stands to 
the left and so on. 



The Proposition |2.1| implies the following formula for the product of n Lie algebra ele- 



ments Vf 



v 1 ...v n = d ei ...d en A ( E cr Bi,.-,« I <(i) • • • v s[ifs(i) ■ ■ • e s( i) 



m=l m - \ 1=1 seS l 



,(2.16) 

ei=...=e„=0 

here is a set of all injections of I numbers 1, . . . , I into n numbers 1, . . . , n which keeps 
the order. A more explicit version of this formula requires a splitting of n elements v j into 
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m sets, with rij elements in each set, m being an arbitrary integer number: 1 < m < n. 
Denote by (s 1; . . . , s m ) an injection of the union of m sets 1, . . . , into the set 1, . . . , n 
which preserves the order within each set 1, ...,nf si(i) > si(j) if i > j. Consider now a 
symmetrized product 

™llz=il#i)! seSm 

here 

Di = ia a C a a a / 1 m ..,^; 1 (2.18) 

and is the number of indices % for which rii = I. The r.h.s. of eq. (2.16) is equal to the sum 
of all such D taken over all numbers m, all possible splittings and all injections (s±, . . . , s m ). 

Let us apply this presentation to the product A(t s ^) . . . A(t s ^) appearing in the second 
line of eq. (2.15). Suppose that we permute some of A(t s ^). A term D coming from a partic- 
ular injection (s\, . . . , s m ) still remains if this permutation does not change the order within 
the rij-element sets into which the elements A(t s u\) are split by the injection. Therefore we 
can combine the integrals over < A(t s m) < . . . < A(t s ^) < 1 which come with a partic- 
ular term D into one integral over the regions < A(t s ^ ni ^) < . . . < A(t s ^ Si (i^) < 1 for 
1 < % < m. This operation leaves a subgroup of S which permutes the images s«(l), . . . , Si(rii) 
for a given i. Its composition with S l n creates a redundant group S. The sum over its ele- 
ments can be removed by relabelling the integration variables and adding an extra factor n\ 
which cancels the same factor in the denominator of eq. (2.15). It is not hard to see that 
what remains is the mth order term in the expansion of the exponential in eq. (2.13). This 



completes the proof of the Proposition 2.2 



A presentation of the holonomy operator as an exponential enables us to use the Weyl 
character formula for the calculation of its trace. For an element v of the Lie algebra su{2) 

^ = sin(aH) = H f *S e M (2 . 19) 



Combining this equation with eq. (2.13) we conclude that 

Tr a Pexp ( I A,dxA = ^# = J*L / (2.20) 
\Jk. J sin \v\ sin \v\ J\d\=§ An \a\ 
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^ = EQ,.„ / dt 1 ,...dt n {A*(t 1 ),...,A-{t n )}. (2.21) 

n=l J0 

This is the exponential presentation of the trace of holonomy that we were looking for. 

We can apply the formula (2.20) to the holonomies along the link components of eq. (1.1): 

n / K d 2 a- 

" ' ' ' ' r " r ' " * (2.22) 









sin 





x exp 



e ^p / d * x A . Q vAp + -A, ■ (A u xA p ))+iKY^ aj 

71 Jm \Z 6 / - =1 



here 



«i = E <£,...,«» / 4 j) • • • dtifW). ■ ■ ■ , ^)}, (2-23) 

n=l J0 

and is a parametrization of a link component £j. We put K instead of A; as a factor in 
front of the integral in the exponent of eq. (2.22) in order to be able to ignore the 1-loop 
corrections to the propagator of the original Chern-Simons theory (1.1) (see e.g. [7] and 
references therein). The formula (2.22) indicates that the terms Sj ■ vj may be considered as 
extra vertices in the quantum Chern-Simons theory. In other words, Feynman rules for the 
quantum theory (2.22) include a propagator (i.e., a Green's function) 



ITT 



(A^ Xl )A b u (x 2 )) = --5 ab n, l/ (x 1 ,x 2 ), 



a usual cubic vertex 



K 



JM 



3ir Jm 



and extra vertices coming from the holonomies along the link components 



V® = iKa 3 ■ C ai _ an 



dt?...dt^{A(t?),...,A(t^)}. 



In particular, 





= iKdj ■ ^ 




= —iKdj 




2 

= ^Kaj ■ 



dt ij) A{t 



dt?dt?{A(t?) x A(t?)}j , 

J^dt?dt?dt? ({(AA J) ) ? At?)) ? 

? ? At?))}) 
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(2.24) 

(2.25) 

(2.26) 

(2.27) 
(2.28) 
(2.29) 



here 

a 3 ■ {A(t?) >< A(t?)} = e abc a^{A b (t?), A c (t^)}. (2.30) 

A symmetric bilocal (l,l)-form (i.e., a 1-form in both variables x, y) Q^ v (x,y) of eq. (2.24) 
is a Green's function of the operator e >wp d v (i.e., of a differential d). It should satisfy an 
equation 

d y n(x,y) = 5 {3 \y-x) + d x n(x,y), (2.31) 

here S^(y — x) is a 3-form ^-function while £l(x,y) is a (0,2)-form, i.e. a 2-form in y and a 
0-form in x. If M — S 3 and S* 3 is presented as 1R 3 with an infinite point, then 

I yP _ x p 

tt^(x,y) = — e^p- jg. (2.32) 

47r \y — x\ 6 

For more information on Q(x,y) see, e.g. [7] and references therein. 

There are some extra vertices coming from the expansion of the prefactors J^i -| into the 
powers of vj. These vertices do not assemble into an exponential. According to the standard 
combinatorics of Feynman diagrams, the trivial connection contribution to the path integral 
over [DAp] in eq. (2.22) can be presented as a product of two factors: 



Z£L (M, C, k) = (1 + G(d u ...,d n ,K)) exp 



J2K 1 - l L^(d 1 ,...,d n ) 



1=0 



(2.33) 



here G(di, . . . , a n ; K) is a sum of all Feynman diagrams containing at least one vertex coming 
from si ^|l.| and L^(di, . . . ,a n ) is a sum of all connected /-loop Feynman diagrams which 
contain only the vertices (2.25) and (2.26). Each of these vertices carries a factor K while 
the propagator (2.24) is of order K' 1 . As a result, /-loop diagrams have a factor K l ~ l which 
we made explicit in the exponent of eq. (2.33). The vertices coming from si |^l| do not carry 
the factor K, therefore the diagrams that contain such vertices have only zero or negative 
powers of K. Thus, combining the Taylor series expansions 

00 

L (1) (ai, . . . , a n ) = L m(d u . . . , On), (2.34) 

m=2 
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[1 +G(di,...,a n ;K)]exp 



1=1 



(2.35) 



1+ J2 K- m P myl {a^ ...,a n ), 



l,m=0 

±j m [CLi y . . . , CL n ) and Pi >m ,(ai, . . . ,a n ) being invariant homogeneous polynomials of order m, 
with eqs. (2.33) and (2.22) we arrive at Reshetikhin's formula (1.8). 

A quadratic polynomial L<i comes from the Feynman diagram containing only one prop- 
agator (2.24) both endpoints of which are attached to link components. Since 

dxi <t> dxjil(xi, Xj) = kj, (2.36) 

lij being the gaussian linking number, we conclude that 

n 

L 2 (ai, . . . , On) = ^2 l ij ®i ' ®r (2.37) 



This completes the proof of the Proposition 1.2 



We will also need in the future the polynomials 



L 3 (a 1} . . . , On) = Hjl ' K' X 5 *)> ( 2 - 38 ) 



i,j,k=l 



it 



L 4 (ai, . . . , a n ) = Hj,ki x %) ' x a/)- (2.39) 

i,j,k,l=l 

It will become clear from the study of Feynman diagrams in Sections |5] and |] why we use 
these particular group weight structures. 

Consider now the group weights associated to connected tree level diagrams contributing 
to the polynomials L m { CL\ , . . . , CL n ). These diagrams are combinations of tree diagrams of the 
original Chern-Simons theory (1.2) sewn by the new vertices (2.26) which are themselves 
tree diagrams (see Figs. 1-4). In both types of tree diagrams the segments are 5-symbols S a b 
and the elementary cubic vertices are proportional to the group structure constants e a bc- 

We may associate an invariant monomial F m (b\, . . . ,b m ) which is linear in all vectors 
bj, I < j < m to & combined group weight coming from a Feynman diagram with m 
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vertices (2.26) by placing vectors bj at the bottom of the tree diagrams associated with 
its vertices (2.26). Then in order to get the actual group weight of the Feynman diagram 
that would contribute to the polynomial L m (di, . . . , a n ), we should substitute n vectors dj 
for m vectors bj depending on which vertex (2.26) comes from which link component Cj. 
Since every tree diagram with more than three external legs contains at least two F-shaped 
configurations with two vectors b attached to each of them (see Fig. 5) and since the cubic 
vertices produce antisymmetric tensors e a b c , we conclude that a polynomial F m (b\, . . . , b m ), 
m > 4 is zero if at least m — 1 of m vectors bj are parallel. It is obvious that the same is 
true for m — 3. 

Proposition 2.3 The polynomials L m (di, . . . , a n ) are produced from invariant homogeneous 
polynomials F m (bi, . . . , b m ) of order m by substituting n vectors dj in place of m vectors bj. 
The polynomials F m (pi, . . . ,b m ), m > 3 are equal to zero if at least m — 1 of m vectors bj 
are parallel. 

This proposition will play an important role in extracting the multivariable Alexander poly- 
nomial from the r.h.s. of eq. (1.8) in [3]. 

Suppose for a moment that the link £ has only one component, i.e. it is in fact a knot K,. 
An immediate consequence of the Proposition |2.3j is that only a quadratic term survives in 
the exponent of the Reshetikhin's formula (1.8). As a result, the whole formula is reduced 
to eq. (1.5). Thus we produced yet another proof of the first part of the Melvin-Morton 
conjecture 1.1. 

3 The Exponent of Reshetikhin's Formula 

In the remainder of this paper we are going to study more closely the structure of the 
exponent of Reshetikhin's formula (1.8). 

Conjecture 3.1 If Li( ) = for all I < m, then the coefficients of the polynomial 

L m (d\, . . . , a n ) are proportional to the mth order Milnor's linking numbers im of the link 
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C: 

L m (a u ...,a n ) = ^ J2 4ri.,i m Tr (^-4)---(^-«i m ), (3-1) 

l<ii,...,i ra <n 

/iere a = (<7i, o"2, 03) is a 3- dimensional vector formed by Pauli matrices. 

We use the scalar product ■ Oj instead of simply in order to stress that we are dealing 
with the su(2) Lie algebra element in the fundamental representation. Note that iff = Z^, 
so eq. (3.1) for m = 2 is consistent with eq. (1.10). 

It is known that Milnor's link invariants of the higher order are not well defined (at least, 
in Z) if lower order invariants are non-zero. The same property is shared by the polynomials 
L m . Their coefficients can not be restored unambiguously from the value of the partition 
function Zg*{ an (M, £, k) because there are changes in integration variables a 3 - which keep 
the form of eq. (1.8) but still alter the coefficients of higher order polynomials L m if the 
lower order polynomials are non-zero. Suppose for example that we substitute a vector 3j 
by a vector a'j obtained by rotating Sj around another vector a k by an angle (pj t k- 

CO ^pTTl 

4 = J2 Pfc x [a* x • • • [o* x • • •] • ( 3 - 2 ) 

m times 

A substitution of eq. (3.2) into the quadratic polynomial L 2 of eq. (2.37) generates, among 
others, new cubic terms, so that 

4? = l U + x ij¥j,k, 1 < i < n, i^j,k. (3.3) 

To put it differently, we need to know the coefficients l\^ k only up to these transformations. 
In fact, as we will see in Section [5|, the Feynman rules of Section ^| predict these coefficients 
only up to this transformation due to the dependence of the vertices (2.26) on the choice 

of a zero point in t^' parametrization of the link component Cj. However this ambiguity 

(3) 

disappears if the gaussian linking numbers kj are equal to zero. Then the coefficients l\j k are 
well defined and turn out to be proportional to the triple Milnor's invariants. 

There are many other possible rotations of the integration variables 3j. For example, 
a vector 3j can be rotated around the vector a« x Sj. This transformation will change the 
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coefficient by an amount proportional to Zy. It will also cause a change in the coefficient 
Pij of the preexponential polynomial 



0,2[ a l; 



i,j=l 



(3.4) 



due to a nontrivial jacobian factor. 

Let us now consider some circumstantial evidence in support of the Conjecture |3.1| . 
Milnor's linking numbers are invariant under tying a small knot on a link component. Let 
L be an n-component knot in a RHS M and /C be a knot in S 3 . We cut an infinitely small 
3-dimensional ball B 3 , whose center belongs to C\, out of M and we cut another infinitely 
small ball B 3 , whose center belongs to /C, out of S 3 . Then we glue the boundaries d(M\B 3 ) 
and d(S 3 \ B 3 ), thus producing a new link £' in M. In other words, we "tie" a small knot 
K, on the component C\ of C Milnor's invariants of £ and £' coincide. 

Let us see what happens to the exponent of eq. (1.8). According to [1], 

l~K 1 



Zgl^(M, £';k) = z£l >an (M,£-,k)Z ai (S 3 ,)C-k^ 



(3.5) 



2 sin (f^ai) 

Combining the formula (1.8) for Zj&\ _ a (M, C, k) with the formula (1.5) for Z ai (S 3 , K\ k) 
we can easily derive the formula (1.8) for Z^*\ a {M, £'j k): 



zS..., an (M,C;k) 

= Z {tr \M; k) 



(3.6) 



" (Kd 2 d- 



exp 



mK 



va 



1+ £ K- m P m ^d Xl ... 



l,m=0 



m=2 



e 2K 



sin (f ) 


7T <2i 






sin(7r|ai|) 



J(K|ai|,K). 



As we see, the exponent remains the same except for the trivial change of framing of L\. 
I'll = hi + v - This provides a confirmation for our conjecture that the coefficients of the 
polynomials L TO ( ) are proportional to Milnor's linking numbers. 
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4 Flat Connections in the Link Complement 



The strongest evidence in support of the Conjecture is provided by the relation be- 
tween the large K asymptotics of the Jones polynomial of a link and the flat connections 
in the link complement. In this section we will assume for simplicity that our manifold 
M is a 3-dimensional sphere S 3 . Consider a large K limit of the Jones polynomial (1.1) 
when the ratios (1.9) are kept constant. Then (see [2] and references therein) the invari- 
ant Z au _ j(ln (S, C; k) can be expressed as a path integral over the connections in the link 
complement: 

Z n „„(S 3 ,£;fc) = JV>AJexp(jrS'cs{AA), t 4 - 1 ) 

S-M = ±TW d 3 x (AAA, - \aAa) (4.2) 

1 n 



Y Tr 
2 ^ 

z i=i 



AAx* 1 N> , A, Ax 



q{3) >* J \Jc { 2 j) 



In this formula and C% are t wo basic cycles on the boundary of the tubular neighbor- 
hood Tub(£j). The cycle C[ is a meridian of Cj, it can be contracted through Tub (£.,). 
A cycle Cf is a parallel, it has a unit intersection number with C[ and it is defined only 
modulo C\ . A self-linking number Ijj is by definition a linking number between Cj and 
■ The path integral (4.1) goes over the gauge equivalence classes of connections A^ which 
satisfy the boundary conditions 

Pexp ^J> (j) A^dx^ = exp (^ a 3 a j^ = exp(i7ra ■ dj) up to a conjugation. (4.3) 

The large K limit of the path integral (4.1) can be found with the help of the stationary phase 
approximation. The invariant Z ai ^ an (S 3 , C; k) will be presented in the form (1.4), however 
this time the sum will go over the flat connections in the link complement S 3 \ J2]=i Tub(£j) 
satisfying the boundary conditions (4.3). On the other hand, the same large K asymptotics 
of Z^ an (S 3 , C; k) can be found by applying the stationary phase approximation to the 
finite dimensional integral over the vectors dj in Reshetikhin's formula (1.8). The invariant 
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will be presented as a sum over the conditional stationary points of the phase 



L m {Si, . . . , 




(4.4) 



m=l 



the conditions being eqs. (1.9). Therefore we conjecture that there is a one-to-one corre- 
spondence between the flat connections in the link complement, which are close to the trivial 
connection, and the stationary phase points of the phase (4.4), so that their contributions 



Let us make this relation more precise. Consider a 1-parametric family of flat connections 
A^(x, r), r > in the link complement, which starts at the trivial connection: A^(x, 0) = 0. 
The holonomies of a flat connection define a homomorphism from the fundamental group of 
the link iri(C) into SU{2). In Wirtinger's presentation of tti(C) the link is projected onto a 
plane and the fundamental group is generated by the meridians Cjj of the pieces into which 
the link components Cj are split by over crossings. If two pieces Cjj and of a link 

component Cj are joined at the overcrossing of Cj by a piece Ck,i of Ck, then the elements 
Cj ti ,Cj ji+ i,Ckj € tti(£) satisfy the equation 



the signs in the exponents depend on the signature of the overcrossing. These relations imply 
that for a given j the holonomies of A^(x,r) along the meridians Cj ti are all equal to the 
leading order in r: 



to Z£l an (S 3 , C\ k) are equal. 




±i 



(4.5) 



k,l ^j,i U k,l ? 




(4.6) 



here a'j are the vectors indicating the directions in the Lie algebra su(2) in which the trivial 
connection of r = in deformed as r grows. 



Consider a knot /C in the link complement S 3 \jyj =1 Tub(£j), we denote as l 0j the linking 
numbers of /C and Cj. It is not hard to see that to the leading order in r 




(4.7) 
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so that if we attach a /^-dimensional representations to /C, then 



Tr^Pexp (J A^x, r)dx^j = (3 



u'=i 



+ 0(r 3 ). 



(4.8) 



Consider now a family of conditional stationary points of the phase (4.4): 

a { ; t) (T) = a"T + 0(T 2 ). 



(4.9) 



Let us add the knot /C as a Oth component to the link C (i.e., multiply the r.h.s. of eq. (1.1) 
by Tr^Pexp {§ K A^dx^)) and see what happens to the contribution of the stationary phase 
points (4.9). The new exponent of the formula (1.8) should include the terms containing the 
vector b corresponding to /C: |6| = (3/K. We assume that (3 ~ 1 as K — > oo. Then, to the 
leading order in r, K and (3, we should account only for the bilinear term b ■ (SjLi hjO-j) 
in the new exponent. As a result, the contribution of the stationary phase point (4.9) is 
multiplied by the factor 

2" 



K d 2 b 
|5|=£ 4vr |6| 



exp 



iirrKb ■ I ^ l^jd"- 



(3 



+ C(r 3 ). (4.10) 



This factor should be interpreted as the trace of the holonomy of the flat connection corre- 
sponding to the stationary point (4.9). Comparing eqs. (4.10) and (4.8) we conclude that 



a! j = d" 3 



(4.11) 



for the family of the flat connections A^(x, r) that corresponds to the family of the stationary 
phase points (4.9). 

Milnor's linking numbers 4fi.,j m allow us to formulate the necessary conditions that the 
vectors a'j have to satisfy so that the flat connections A^{x,t) with the holonomies (4.6) 
exist. Let us briefly review the algebraic definition of these numbers. Consider a parallel Cj, 
by definition it is an element of ni(£) which is homologically equivalent to the parallel 
and also commutes with the meridian C^\. The element Cj can be expressed as a product of 
powers of the meridians Ck,i since they generate the whole group ni(C). Milnor showed [8] 
that the parallel Cj can be expressed only in terms of the meridians Ck,i (one meridian 
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per link component) modulo the elements of the lower central subgroup n[ q \C) of 7Ti(£) 
(^(C) = 7Ti(£), ir[ n+1 \C) = [n[ n \ tvi(C)]) for any arbitrarily big value of q G Z. Suppose 
that the first non-zero Milnor's numbers appear at order m. Then we choose q > m. It is 
easy to see that for the family of flat connections (4.6) the holonomy along the elements of 
n[ q) is equal to 1 up to the order r q . We are interested in the holonomy along Cj up the 
order r m_1 , so we can neglect the elements of n[ 9 \£) and use the expression for Cj in terms 
of Ck,i modulo tt[ 9 \jC). Then, according to eqs. (4.6) and Milnor's definition of l^l i m as the 
coefficients of the Magnus expansion of this expression, 

Pexp (j>_ A^x,r)dx^ (4.12) 

= i + h- 1 £ i&^A* <)•••(*• «L-J + o(r m ) 

l<ii,...,i m -i<n 

(we actually used the relation l iti = exp[27rr<? ■ a*] coming from eq. (4.6) rather than d t i = 
1 + Xi, Xj being an indeterminate, which is a standard form of the Magnus expansion, see 
also [10]). Since the parallel Cj commutes with the meridian Cj t i, we come to the following 
conclusion: 

Proposition 4.1 If the 1 -parametric family of flat connections A ^(x, r) defined by eq. (4-6) 
exists in the link complement, then the vectors a'j satisfy the condition 



l<ii,...,j m _i<n 



0, l<j<n. (4.13) 



Let us compare eq. (4.13) with the stationary phase condition for the integral of eq. (1.8). 
Suppose that Lj( ) — for I < m. Then eq. (4.9) presents a family of conditional 

stationary points of the phase (4.4) if 

dL m (a", . . . , a") „ . . 
da] X a > = °" (414) 

Then it follows from the invariance of Milnor's linking numbers im under a cyclic per- 
mutation of the indices that 

Proposition 4.2 The identification (4-H) of the flat connections (4-6) with conditional 
stationary points (4-9) of the phase (4-4) ^ s consistent with the conjectured expression (3.1) 
for the polynomials L m (ai, . . . , a n ). 
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This proposition supports the Conjecture |3.1| but it does not help us to fix the coefficient 
in front the trace in eq. (3.1). This can be done by comparing the dominant part of the 
stationary phase of eq. (1.8) 



-T m L m (a' 1} . . . , a' n ) 



(4.15) 



with the Chern-Simons action (4.2) of the flat connection A^{x,r). The action S' cs has the 
following property (see, e.g. [9]): if the connections and A^ + SA^ are both flat, then 



3 , ca [A ll + 5A lt ]-S , cg [A ll \ = - y £Tt 



, s SA u dx^ 

b) ^ 



c: 



A^dx" 



(this is a general property of a classical action: 5 {j^\pq — H(p, 
other hand, eq. (4.12) implies that 



p5q 



*2, 



(4.16) 



. On the 



c. 



b) 



A^Xjrfdx 11 = (wrr) 



m— 1 



E 



l<il,...,i m _i<n 

Combining eqs. (4.16) and (4.17) we conclude that 



dr 



t* e e 



(4.18) 



l<ii,...,i m <n 

After integrating this equation over r we arrive at eq. (3.1). 



5 A Triple Milnor's Linking Number 

Milnor's invariants ijf ... i m can be expressed as integrals of differential forms constructed with 
the help of the Massey product (see [13], [14] and references therein, a simple introduction 
into this subject together with the relevant formulas can be found in [15]). We are going to 
check the Conjecture [TT] for the polynomials L3 and L4 by comparing the Feynman diagram 
formulas for their coefficients with these expressions. 

Preliminaries 

We start by introducing some useful notations. Let £ be an n-component link in a 
RHS M. Suppose that we cut out a tubular neighborhood Tub(£j) from the manifold M 
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and then glue it back after switching its meridian C[ j ^ and its parallel C 2 \ We call such 
procedure an S-surgery on Cj. With a slight abuse of notations, we denote as Tub' (£_,•) 
the tubular neighborhood when it is glued back as a result of S'-surgery. We denote 
as ^i 1 ...i m ,j 1 ...j i the manifold constructed from M by removing the tubular neighborhoods 
Tub(£j 1 ), . . . , Tub(£j m ) (we will also assume these neighborhoods to be infinitely thin) and 
performing S'-surgeries on the link components C^, . . . Cj r 

We denote as Uj a closed 1-form defined in Mj by a condition 
This form can be expressed with the help of the Green's function (2.24): 

Wj -(.) = jf j n(*w (5-2) 

here t^) is a parametrization of Cj and we slightly abused the notations by using instead 
of x(t^ as the argument of f2. The linking numbers kj can be expressed with the help of 
the forms Uj. 

Uj= I Uj= I Ui= I I n(t {i) ,t {j) ). (5.3) 

Another useful property of Uj is that if u is a smooth 1-form in an infinitely thin tubular 
neighborhood Tub(Cj), then 

/ Uj A u = <p uj. (5.4) 

JdTnb(Cj) JCj 

The following object appears naturally in the formulas for Milnor's invariants. Let u±, U2 
be two 1-forms defined on a knot K, parametrized by < t < 1. An "iterated commutator" 
[^1,^2] is a bilocal (l,l)-form 

[u l ,u 2 ]{t 1 ,t 2 ) = s\gn{t 1 -t 2 )u 1 (t 1 )u 2 (t 2 ). (5.5) 

If both forms U\ and u 2 are multilocal, then in our notations 

„ „ m+n 

/ [<^i, ^2] = / TJ dtj u 1 (t 1 , . . . , t m )u 2 (t m+ i, . . . , t m+n ) 

J min(£i , ...,t m )>max(t m +i , ...,t m + n ) j—\ 

„ m+n 

— J F[ dtj ui(ti, . . . , t, m )u 2 {t m+ i, . . . , t m+n ). (5.6) 

Jmin(t m+ i,...,t m+n )>max(ti,...,t m ) - =1 
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Obviously, the definition of the iterated commutator (5.5) depends on the choice of the 
zero-point of t parametrization. If this zero-point is shifted by At, then 

lim ^Mp^A = 2^(0) f u 2 - 2u 2 (0) I Wl . (5.7) 
At^O At Jk Jk 

Therefore if both integrals § K uji j2 are equal to zero, then the integral Jici^i, ^2] is well defined. 
The iterated commutator appears in our calculations due to the following 

Proposition 5.1 Let JC be a knot in a manifold M, C\ y2 being the meridian and parallel on 
<9Tub(/C). Let M' be a manifold constructed by S-surgery on JC. Let u be a closed form in 
M \ Tub(/C) satisfying a condition 

I uj = 1. (5.8) 

Let uj±,uj 2 be two closed 1-forms in M and suppose that 

j> uo = j> uj\ = j> uo 2 = 0. (5.9) 

Then the forms uji,uj 2 and uj can be extended into M' . If the tubular neighborhood Tub(/C) 
is infinitely thin then 

[ uj 1 A uj 2 A uj = i / [uji,uj 2 ]. (5.10) 
■/Tub'(/C) 2 Jk 

To prove this proposition we introduce the functions f\ and f[ such that 

df x = ut inside Tub(/C), (5.11) 
df[ = u x inside Tub'(/C) (5.12) 

and they coincide on the common boundary: 

/ilaTub(/c) = /llaTub'(^) • (5.13) 



Then 



/ lui A u) 2 A uj = / df[Auj 2 AuJ = f\uj 2 A uj = — f f\uo 2 (5.14) 

./Tub'(/C) " JTub'(/C) " JdTuh(K) JK. 
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and 



[ui,u 2 ]= dt 2 oj 2 {t 2 ) / dtiui(ti) - / dtiutfa) = -2 f 

JK JO Ut 2 JO J JK 



-2 f ftu 2 . 

K. 



(5.15) 



This proves the proposition. 
Milnor's Invariant 

Suppose that the following linking numbers are all equal to zero: 

hi = 0, 1 < i, j < 3 



(5.16) 



(a condition Ijj = 0, 1 < j < 3 is not necessary to define the triple Milnor's invariants but 
it will simplify our formulas). These conditions allow us to extend the forms ui,uj 2 , uj 3 into 
the manifold Mjjg. The triple Milnor's linking number of the link components C\, C 2 , £3 is 
equal to the intersection number 



123 = / ui Au 2 Au 3 . 

A% 3 



(5.17) 



The r.h.s. of this formula can be expressed as an integral in M. Indeed, we may split the 
integral (5.17): 
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^123 = / UJi A LU 2 A LU 3 + ^ UJi AUJ 2 A UJ3. 

JM 123 ^ ./Tub' (Cj) 



(5.18) 



If the tubular neighborhoods Tub (£7) are infinitely thin then the first integral in the r.h.s. 
of eq. (5.18) is equal to 



uj\ A lo 2 A oo 3 = Y123 = / A dJ 2 A CU3, 
M123 Jm 



(5.19) 



while the integrals in the sum can be transformed with the help of the Proposition |5.1| : 

f u>i A ujj A u k = X^ k = - [ [uji, ujj], i^J^k. (5.20) 



Therefore 



i(p) 

t 123 



1 123 "r ^12,3 "T ^31,2 "I" ^23,1 



'(7) 



(V) 



'(7) 



(5.21) 
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Feynman Diagrams 

Now we will use the Feynman rules derived in Section |2| in order to calculate the cubic 
coefficient li 23 . The Feynman diagram contributions to ^23 are depicted in Figs. 6 and 7 
up to the permutations. The diagram of Fig. 6 contains three propagators (2.24), one cubic 
vertex (2.25) and tree vertices (2.27). Its contribution to the exponent of eq. (1.8) is equal 
to 



r,(6) 
-^123 



-2in 2 Ka 1 ■ (a 2 x a 3 ) / d 3 y e UlU2V ' A (5.22) 

JM 



xf dxf n fllUl (x 1 ,y) f dxf Vt^ 2V2 {x 2 ,y) f dx^VL^ V3 (x 3 ,y), 
or, in view of eq. (5.2), 

D{% = -2m 2 Kax ■ (a 2 x d 3 )Y^. (5.23) 

A diagram of Fig. 7 consists of two propagators (2.24), two vertices (2.27) and one ver- 
tex (2.28). Its contribution is equal to 

■°8i = -i^Kai ■ (a 2 x a 3 ) / f dx 2 Q(x 2 ,-), f dx 3 tt(x 3 ,-) , (5.24) 

or, in view of eq. (5.2) 

£>8i = -™ 2 Kd x ■ (a 2 x d 3 )xg\. (5.25) 
A corresponding cubic term in the exponent of eq. (1.8) is 

ZmKlf^ax ■ (a 2 x a 3 ) (5.26) 

so that 

^123 = "g 71 " (^123 + -^12^3 + ^31,2 + X 23 ^j . (5.27) 

Comparing this expression with the formula (5.21) for the triple linking number we see that 
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Proposition 5.2 If the linking numbers lu,ljj,lkkikj,ljk,kk are a H equal to zero, then the 
cubic coefficient l\^ k in the exponent of Reshetikhin's formula (1.8) is proportional to the 
triple Milnor's linking number 



The obvious symmetry 



hjk ~ Ijik (5.29) 



leads to the following relation: 



lijl • Si) (a ■ Sj)(a ■ S k ) = 2i l ijk x ' °>k- (5-30) 

After substituting eq. (5.28) into eq. (3.1) and using eq. (5.30) we recover eq. (2.38). This 
proves the Conjecture |3J] for m = 3. 

Note that the contribution (5.25) of the Feynman diagram of Fig. 7 is proportional to 

(7) 

the integral ^ of eq. (5.20) whose value depends on the choice of the zero-point it the 
parametrization t\ of the link component L\. However it follows from eq. (5.7) that this 
ambiguity can be compensated by the change of integration variables (3.2). 



6 A Quartic Milnor's Linking Number 

Preliminaries 

We begin with establishing another formula for the triple Milnor's linking number by 
splitting the integral (5.17) in a different way: 

Hil — u i A u i AWjjl / u>i A u)~ A uj k . (6.1) 

All intersection numbers of the closed 2-form ui A uij in the manifold Mjj are equal to zero, 
so it is exact: 



Ui A ujj = duJij . 
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(6.2) 



Therefore 



/ uoi A uOj A uok = I duoij AuJk = — A (Vh — V Wa (6.3) 

and 



4 V' J ' 2 1 



= /. + i [<*><",■]). ( 6 - 4 ) 



If we introduce the functions /j^, ■ satisfying equations 

dfij = Ui inside Tab(Cj), (6.5) 
df-j = Ui inside Tub'(£j)» (6.6) 

then the integrals of the iterated commutators can be presented in the way similar to 
eq. (5.15): 

/ [ui,Uj] = -2 <b f itk u j = 2j f jt kfVi= f (fj,kVi - fi,kVj) ■ (6-8) 

A 1-form Uij is defined by eq. (6.2) in the manifold Mjj only up to an addition of a closed 
form (that is either ui or ujj). One possible candidate for can be obtained with the help 
of the propagator (l,l)-form Q(x,y): 

1 r r „ 1 



UiM= d 3 xn(x,y) AWj(x) AUj(x)-^ I M;y),Wj(-)] + ^ I M-,y),u t (-)] . (6.9) 

■'Mi, 2 C, I JCj 

Indeed, it follows from eqs. (2.31), (5.4) and (6.8) that 

duij(y) = Ui(y) A uj(y) - / d 2 x Cl(x,y)ui(x) A Uj(x) (6.10) 

JdMij 

-\ fjdn^y)^.)} - £ di« 5(y - s(fM)) / o <W 



+-/ [dn(-,y),a; i (0]+/ dt® 6(y - z(t®)) <ft<^(x(f?>)) 
z Jc 3 Jo Jo 

= Ui (y) A Uj (y) - £ dt® 6(y - x(t®)) dtfcuj(x(tf)) 



+ / dt® 5(y - x(t®)) / dtjpuMtP)) 
Jo Jo 
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— / d 2 x Q(x, y)uji(x) A uoAx) — / d 2 x y)tUi(x) A uoAx) 

JdTub(Ci) JdTub(Cj) 

+ / f2(-,j/H-(-)- / f2(-,J/)wi(0 
= Ui (y) A Uj (y) - £ cftW % - *(*«)) j^" dtfu^f)) 

/ dt<» 6(y - x(t<»)) / (ftpW^*? ))- 

JO JO 



+ 

The last two terms of this equation are equal to zero when y e M^-, however they are useful 
in deriving the formula for contour integrals of Uij by using Stoke's theorem: 

Proposition 6.1 Let C be a contour in which is a boundary of a 2- dimensional ori- 
entable surface E e M. Let {Pj} and {Pj} be sets of points where X intersects Ci and 
Cj, sign (Pj) and sign (Pj) being the signatures of these intersections. Then for ujij given by 
eq. (6.9) 

j uJij = / sign (Pi) / Uj + J2 si S n ( P j) / w »- ( 6 - n ) 

J C J S p JO p JO 

Let us choose an infinitely small meridian c[ l \t^) near the point t^ as the contour C 
and a small disk which intersects Li at the point as the surface X. Then the first integral 
of eq. (6.11) is infinitely small because the singularity of the form u>i near £j is compensated 
by the jacobian measure factor in polar coordinates. 

Corollary 6.1 If the meridians Ci\t^) and C[ j \t^) are infinitely small, then 
h^(td))" lJ Jo ^' Jc[ j) (tU)) 



ft 1- ' r rt u ' 

"«=-y ^ u im ^ = -L Ui - (6 - 12) 



Milnor's Invariant 

Suppose that in addition to the gaussian linking numbers, the triple Milnor's invariants 
are also equal to zero: 

= 1^1 = 0, l<i,j,k<4. (6.13) 

Then the 2-forms UJi A cuj are exact in the manifold Mjgga because the intersection num- 
bers (5.17) are equal to zero. Therefore we can introduce the 1-forms ujij of eq. (6.2) in 

26 



Mi23i- Their values in M 12 34 G Mjggj may be given by eq. (6.9). At this point we may use 
the Massey product (see e.g. [12], [15] and references therein). The following 2-form 

Uij,k = Uij AW3- T^fci A IQj - A io { (6.14) 

is closed and the integral 

4\2 = / Ui ^ k A Ul ( 6 - 15 ) 

is a link invariant which is related to the quartic Milnor's linking number: 

lijll = g (lijjl ~ ljk,il)- (6.16) 
An equivalent but more symmetric formula for l^ k \ can be obtained by integrating by parts: 

AM) 



M = / (6.17) 



K 3 

' 1 1234 
1 



id. 



,ki = o AWfcAwii w fc/ AWiA (6.18) 



2 

1 , , 1 



--j{u k i A A uji + u;^ A ^fe A a;,) - -(wjfc A cjj A cj,- + lu u A A u k ) 



AM) 
ij,kl 

AM) _ AM) _ AM) _ AM) 



We will work directly with the invariants l\j kl because they have the obvious symmetries: 



L ij,kl — L ji,kl — L ij,lk — l kl,ij- [0.1V) 

They also satisfy a Jacobi identity: 

AM) AM) AM) _ „ , fi 
f 12,34 + f 31,24 + '23,14 ~ u - V°- ZU J 

which indicates that the space of quartic invariants of a 4-component link is 2-dimensional. 

Our goal is to express the invariants l^ k \ as integrals in M. We will work with a particular 
invariant ^^34 in order to simplify our notations. We split the integral (6.17): 

/ ^12,34 = / ^12,34 + J2 ^12,34- (6.21) 

Jm^u JM 12 -ii j=1 JTub'iCj) 

The integral Jm 1234 ^12,34 becomes j M ^12,34 in the limit of infinitely thin tubular neighbor- 
hoods Tub(£j). For any 3-form A u k A u>i we can use eq. (6.9) in order to obtain the 
formula 

j M Wy A u k A W , = + Xg2, - Xjjji, (6.22) 
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here 



Y: 



(7) 



X 



ij,kl 



(10) 



d A y[ I d 6 xuji(x) AoJj(x) Att(x,y)\ A oo k (y) Aui(y), 

M \JM 



i,j,kl 2 



M 



d 3 y[ / [n(;y),Ui(-)] \ Au k (y) Aui(y). 



(6.23) 
(6.24) 



As a result, 



f . . _ 7 (o) 1 7 (o) 1 7 

/ ^12,34 — -^12,34 ~~ 0^31,24 ~~ ^ 



(0) 

23,14' 



(6.25) 



here 



7 (0) _ V (7) , 1 /^(10) y(10) y (10) y< 

^y.Jfel ~~ J y,A;Z ~r r> l/Hj,^ ^i.i.JW ~r k,l,ij 



(10) 



(6.26) 



Next we turn to the integrals over Tub'(£j) in eq. (6.21). Consider, for example, 
/Tub'(£ 4 ) ^12,34- After integrating by parts, we can turn this integral onto a sum 



/ ^12,34 — h + I2 

J Tub' {d) 



h = f (vi2 Au 3 - -u 31 Au 2 - -uj 23 A wi) A o; 4 , 

JTub'(£ 4 ) V 2 2 / 

J 2 = - / fo;i2 A u; 34 - -u; 3 i A u; 24 - -o; 23 A w 14 

2 JaTub(£ 4 ) V 2 2 



(6.27) 
(6.28) 
(6.29) 



To calculate I\ we rewrite it as 

•/Tub'(£ 4 ) LV 2 2 / 



'(A) 

2(^31 + /l,4^3) A ^ 2 - -(^23 - /2,4 W 3) A Wi 



(6.30) 



A 



All three 1-forms 



^12 — 2^,4^2 + ^2,4^1, ^31 + /l,4^3j ^23 ~ 72,4^3 



(6.31) 



are closed. They also satisfy the condition (5.9) of the Proposition [54] for fC = £4 because 
of eqs. (6.4), (6.8) and (6.13). Therefore we can apply the Proposition |54] to the r.h.s. of 
eq. (6.30): 



Vl2 + ^[^l,^], ^3 



~[uJ 3 l,UJ 2 ) ~ 2^23,^1] 



(6.32) 
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At the same time, Corollary |6.1| allows us to express I2 also as a contour integral. Since the 
form Uij is nonsingular in Tub(£;), then in the limit of infinitely thin tubular neighborhood 



ujh A Uki = - / dt (l) 



Ukl — —-, 



(6.33) 



so that 



l2 = ~\j c ([^12^3] - 7^31,^2] - 2^23,^1]) • (6.34) 
Combining eqs. (6.27), (6.32) and (6.34) we conclude that 

/ ^12,34 = 7 / ( [^12,^3] - ^31,^2] - -[^23,^1]) + 7 / [[^1,^2], ^3], (6-35) 

./Tub' (£4) 4 J£ 4 \ 2 2/4 J£ 4 

or after substituting eq. (6.9) for u>ij, 



f 

JTub' 



Tub'(d) 



7 (0 1 7 (0 _ 1 7 (0 , 3 V-(9) 



(6.36) 



here 



v-(8) 
ij,k,l 



(0 _ 1 /V(10) _|_ v( n ) 

[^(-,*),^(-)]V^W 

[[cjj,^-],^]. 



(11) 



i,j,k,l 



2 
1 

6 ./£, 



(6.37) 
(6.38) 
(6.39) 



It remains now to put together eqs. (6.17), (6.21), (6.25) and (6.36): 



(AO 



1 



~.Zki,jl 



1 



ij,kh 



(6.40) 



here 



,(0) 



+ + ^ 



0') 



ij,kl ' ij,kl ' ij,kl ' ij,kl 



7 W 

J ij,kl 



(6.41) 



v (7) y (10) y (10) y(10) y (10) 

J ii,fei ' ^i,j,kl ^j,i,kl ^k^ij ^l,k,ij 

_^( n ) 4_ ^( U ) _ ^( n ) _i_ 

^j,i,l,k i,j,l,k ^j,i,k,l "r ^-i,j,kp 



Z' 



ij,kl 



y(9) 



v(9) 
ij,l,k 



X, 



(9) 



kl,i,j kl,i,j' 



X 



(9) 



(6.42) 
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Feynman Diagrams 

The Feynman diagrams that contribute to the quartic term 

(Si x Uj) ■ (3 k x 3i) (6.43) 

are drawn in Figs. 8-11 up to permutations. Their contributions are easy to calculate with 
the help of the Feynman rules derived in Section 0: 

D% = Un 3 K(3i x Sj) ■ (3 k x 3 t )Y^ (6.44) 

D® k j = Am 3 K(3i x 3j) ■ (3 k x a,)*gl,i, (6-45) 

= U^K(3i x a,) • (3 k x S,)^Sh' ( 6 ' 46 ) 

= -^K(a t x a,) • (3 k x (6.47) 

In eq. (6.45) we kept only the part of the contribution of the diagram of Fig. 9 which is 
proportional to the quartic term (6.43). The sum of all these diagrams with appropriate 
permutations should be equal to the quartic term in the exponent of Reshetikhin's for- 
mula (1.8): 

4iirKl\f kl (3i x 3j) ■ (3 k x 3i). (6.48) 

Therefore comparing eqs. (6.44)-(6.47) with eqs. (6.41) and (6.42) we conclude that 

tf! kl = n 2 (Z tlM + Z^ kl ). (6.49) 

The coefficients lf^ kl as defined by this equation, have the symmetry properties (6.19), how- 
ever they do not satisfy eq. (6.20). This can be fixed with the help of the Jacobi identity 

(3i x 3j) ■ (3 k x 3i) + (3 k x 3i) ■ (3j x 3i) + (3j x a k ) ■ (a, x 3 t ) = 0, (6.50) 

which allows us to add the same quantity to the coefficients /j^, 41j7 anc ^ ljk,u ( an d the ones 
obtained by permutations (6.19)) without changing the value of the sum (2.39). Choosing 
this quantity to be 

(Zij,ki + Zkiji + Zjk,u) (6.51) 
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we get a new set of coefficients 



C = * 2 



(6.52) 



which satisfy the analog of the Jacobi identity (6.20) 



4/',W + ^li,jk + l\l,ik — (6.53) 

in addition to the symmetries 

Z (4) --/ (4) --/ (4) -/ (4) f6 54) 

l ij,kl ~ L ji,kl — L ij,lk — L kl,ij- {U.O<±) 

Finally, comparing eqs. (6.52) and eq. (6.40) we arrive at the following 

Proposition 6.2 If all the gaussian linking numbers and triple Milnor's invariants involving 

the indices i,j,k,l are equal to zero, then the J^th order coefficients Z- 4 ^ of the exponent of 



Reshetikhin's formula (1.8) normalized by the symmetries (6.54) an d Jacobi identities (6.53) 

l ij,kl' 



are related to the quartic Milnor's linking numbers 



hl,M — o (hj,ki ljk\u) (6.55) 



7T 

The symmetries (6.54) lead to the equation 

E C'/.uu;' 1 ' 1 '^ • •••(*• %J = -2 E !ju,j:,.n (( 'j> x % a ) • x %J(6.56) 

Substituting it together with eq. (6.55) into the r.h.s. of eq. (3.1) we obtain eq. (2.39). Thus 
we verified the Conjecture |3.1| for m = 4. 



7 Discussion 

Reshetikhin's formula (1.8) seems to be an excellent tool for the study of links. It separates 
the exponent, which is of order K 1 , and the preexponential factor, which is at most of order 
K°, for future use in stationary phase approximation. It also incorporates naturally the 
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Feynman diagrams of links and seems to relate them to Milnor's invariants. This result may 
not be so surprising in view of the fact that, as it was demonstrated in [10] and [11], Milnor's 
linking numbers are Vassiliev's invariants of the link. The ambiguities of link Feynman 
diagrams also seem to match the ambiguities of Milnor's linking numbers. Both of these 
ambiguities appear to be related to the freedom of changing the integration variables in 
Reshetikhin's integral (see eqs. (3.2) and (3.3)). It is remarkable that Feynman diagrams 
of the Chern-Simons theory are somehow connected to the Massey product. 

The set of Feynman rules that we derived in Section |2] is not unique because of Lie algebra 
Jacobi identities. In fact, it could possibly be improved in order to incorporate some natural 



combinatorial symmetries. This may facilitate a proof of Conjecture |3.1| about a relation 
between the exponent of Reshetikhin's formula and Milnor's linking numbers. 

The exponent of Reshetikhin's formula might deserve further study. Its stationary phase 
points are in one-to-one correspondence with the flat connections in the link complement 
and the determinant of the second derivatives of the phase is related to their Reidemeister- 
Ray-Singer torsion (it is proportional to the dominant part of the torsion in the limit of 
small phases of the holonomies along the meridians of link components). We use this fact 
in [3] in order to calculate the multivariable Alexander polynomial of the link in the way 
that generalizes the recent Melvin- Morton conjecture [4]. 
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